Importance of the Doppler Effect to the Determination of the Deuteron Binding 

Energy 



G^ 



< 



o 



X 



Yongkyu Ko *, Myung Ki Cheoun and Il-Tong Cheon 
Department of Physics, Yonsei University, Seoul 120-749, Korea 
(February 9, 2008) 

The deuteron binding energy extracted from the reaction ^ H{n,'y)^ H is reviewed with the exact 
relativistic formula, where the initial kinetic energy and the Doppler effect are taken into account. 
We find that the negligible initial kinetic energy of the neutron could cause a significant uncertainty 
which is beyond the errors available up to now. Therefore, we suggest an experiment which should 
include the detailed informations about the initial kinetic energy and the detection angle. It could 
reduce discrepancies among the recently reported values about the deuteron binding energy and pin 
down the uncertainty due to the Doppler broadening of 7 ray. 
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The deuteron binding energy is one of the most important physical quantities in nuclear physics. This plays a 
central role in the determination of the neutron mass, and provides a critical reference energy for the determination of 
high-energy 7 rays measured on the atomic-mass scale 0] . There are several articles to have measured and improved 
the deuteron binding energy more precisely than ever . Some of them have reported their values as precisely as 
up to the order of several electron volts |^,^,^^. Usually, the deuteron binding energy is determined by adding the 
recoil energy of the deuteron to the measured energy of the 7 ray which is emitted from the neutron capture reaction, 
^H(n,j)^H. In the conventional non-relativistic treatment, there seems to be nothing to correct any more, because 
I the initial kinetic energy of the neutron is assumed to be so small that it can be neglected safely. Actually, the kinetic 
' energy of the neutron used in Ref. |Q| is 0.056 eV and that of the thermalized neutron from the neutron source in Ref. 
l^jilHl is assumed to be 0.025 eV at room temperature. 

However, if the two body collision process is treated relativistically, we demonstrate that such a small kinetic energy 
of the neutron causes a considerable uncertainty in the deuteron binding energy because of the Doppler effect. Even 
2 ' the very small neutron kinetic energy gives rise to a significant velocity of the n-p system in the laboratory frame. 
I . As a result, this moving source of the 7 ray can be an origin for the Doppler effect. Naturally this Doppler effect 
Q ' depends on the angle of the 7 ray detector with respect to the velocity of moving source. Therefore, the measured 7 
^ . ray energy is expected to have an unavoidable uncertainty originating from this angle dependence. As will be shown 
^ ' later on, our estimation for the uncertainty shows about 25 eV maximally for the neutron initial kinetic energy of 
J> 0.056 eV and 14.0 eV for that of 0.025 eV, respectively. These values are so large compared to the reported error 2.3 
eV Q, which comes from the 7 ray detector itself. 

Let us suppose that two particles are initially at rest within an appropriate interaction distance, and they come 
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' together to form one body system by emitting a photon. In this situation, the center of momentum frame |T^, which 
is an exact nomenclature in relativistic kinematics but used as the center of mass frame hereafter, coincides with the 
laboratory frame. From the energy and momentum conservation, we get the following relation 

Pi+P2^q + k, (1) 

where p.^'s are four momenta of the initial two particles while q and k are those of the resultant composite particle 
and the photon, respectively. The time component of the momentum accounts for the energy conservation, and its 
space components stand for the momentum conservation. The square of Eq. is Lorentz invariant: 

{pi+P2r^{q + kr. (2) 
In the case of the initial particles at rest, this gives the following equation in the center of mass frame 



(mi -I- ms)^ = + 2cj2 + 2ujy/AP + uj^, (3) 
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where M and rrii are masses of the composite and the initial particles, respectively, and uj is the energy of the emitted 
photon. Solving for w, we obtain the exact equation for the energy of the emitted photon in terms of the masses of 
the initial and the final particles: 



(mi + m2 + M)(mi + r7i2 — M) 

(4) 



2(toi + 7712) 
(Ei + 2M)Eb 



2{Eb + M) ' 

where the binding energy Ei^ is defined as the mass difference between the initial particles and the final product 
Eb = mi + m2 — M . From Eq. we get the binding energy in terms of the produced mass and the energy of the 
emitted photon: 



Eb=uj + yj^"^ + M2 _ M 

where the last expression shows the non-relativistic approximation which is used to determine the binding energy of 
the deuteron so far 

In the case that the incident particle has a kinetic energy and comes to the target at rest, the laboratory frame 
differs from the center of mass frame. Using the Lorentz invariance of the square of four momentum, we obtain the 
following equation 

{V\ + V2fiah = {Pi + P2)i„i = {q + k)i^, (6) 

which is calculated as 



ml + 2m2Eu = Af ^ + 2uj'^ + 2uj^J M"^ + w^. (7) 
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Notice that lo is the energy of the photon emitted in the center of mass frame and is given by 



2-\/m^ + m2 + 2m2Eii 
Here the velocity of the center of mass frame is calculated as 

fJc.m — , -n^ ■ 

m2 + Ell 



ml - AP + 2m2Eu 



(9) 



The above energy of the photon should be transformed to that of the laboratory frame in which the observation is 
carried out. This transformation is nothing but the Doppler effect and is given, in relativistic theory, by 

^.=^/^Hl^, (10) 

where n is the unit vector in the direction of the photon beam. Thus, we obtain the energy of the photon measured 
in the laboratory frame: 



UJl=UJ 



1 - Vc.m COS 

ml + ml- M"^ + 2m2Eii 
2{m2 + Ell - cos6l) 

{Eb + 2M)Eb + 2m2Ki 
2{Eb + M + Ki- |Pi,|cos( 
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where the kinetic energy is defined by Ki = Ei — mi. It should be noted that ui depends on 6, if the neutron beam 
direction is fixed. Alternatively, the same expression can be derived directly from Eq. (|^) in the laboratory frame. 
We consider here three special cases for easy understanding as explained in Ref. |l^. The first case, 9 = which 
represents the transverse Doppler effect, is when we measure the frequency of the photon emitting perpendicular to 
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the incident beam. The second case, 9 = n, which corresponds to the longitudinal Doppler effect in which the source 
is receding, is when we measure the photon emitting anti-parallel to the incident beam. The last case, 9 = 0, which is 
also longitudinal but the source is approaching, is when we measure the photon emitting parallel to the incident beam. 
Solving for the binding energy, we obtain the following equation for the binding energy without any approximation 

Eb^LJi- M + y^M2 + - 2(m2 - uji)Ki - 2uji\p^i\ cos 9, (12) 

where \pii \ = yjKi{Ki + 2mi). This equation, which can be reduced to Eq. (||) as the initial kinetic energy Ki goes 
to zero, contains the initial kinetic energy of the neutrons and the detection angle. Both of them were not considered 
in the non-relativistic analysis of the experiment of ^H{n,j)^H reaction [||. 

In the following, we show the detailed numbers coming from the above 9 dependence and compare to the reported 
error. The measured wavelength of the 7 ray emitted from the n-p capture has been reported as 

Xnp = 5.576 698 8 (55) x IQ-^^m. (13) 

The error in parenthesis comes from the Bragg angle and the lattice-spacing of crystal which are related with the 
equation nX = 2dsni9Bragg- Using the conversion constant fic, we obtain the energy of the photon: 

2-kT%c , , , , 

Lo = — — = 2.223 255 2 (23)MeV. (14) 

The physical constants used in this letter are Q 

he = 197.327 053 (59)MeV fm, 

deuteron mass M = 1875.613 39 (57)MeV, 

neutron mass mi = 939.565 63 (28)MeV, 

proton mass ma = 938.272 31 (28)MeV. (15) 

Since a deuteron mass Af in Eq. (|l|) causes an uncertainty in summing the significant figures, we expanded the 
equation in order to visualize the cancellation of the uncertainty due to the deuteron mass as follows 



ujf {m2-uJi)Ki uJi^jKi{Ki + 2mi) 

Eb = LOi + -— !- - ^ — -— cos 9 

2M M M 

= 2.223 255 2 (23) + 0.001 317 7 - 0.000 000 03 - 0.000 012 2 cos 9 MeV. (16) 

In this result, the second term is the deuteron recoil effect, the third term comes from the initial kinetic energy itself 
and the fourth term stands for the Doppler effect due to the moving source of the center of mass frame. Since the 
authors of Ref. [0 have used the neutron flux of which distribution is approximately Maxwellian with a peak at 1.2 
A(= 0.056 eV), we calculate the binding energies for the three cases of = 0, ^, tt: 

Eb = 2.224 560 1 (23) MeV for cos 6* = 1, 
Eb = 2.224 573 2 (23) MeV for cos 6* = 0, 

Eb = 2.224 585 2 (23) MeV for cos9 = -1. (17) 

This Doppler effect will be reflected as a line width, a line shift of 7 ray or both of them in the experiments. As we 
know from the above numerical calculations, the Doppler effect due to the initial kinetic energy of the neutron causes 
an uncertainty by around 25 eV, which is completely out of the range of the above error 2.3 eV. This error may be 
the best value of the present measurement technique of the 7 ray in that energy scale. Consequently, it is inescapable 
to take the uncertainty due to this Doppler effects into account in the analysis of the ^H{n, ^)'^H experiments. 

In the actual experiment using the thermalized neutron, however, the dependence of the angle is not detected, but 
manifests itself as a line width in 7 ray spectrum because of the following reason. Since the incident neutrons may 
be distributed at random and isotropic in this experiment, we can assume that they are uniformly distributed over 
the whole solid angle An with respect to the direction of the detected photons. Therefore the 7 ray detector fixed at 
the given angle sees all 7 rays triggered by the incoming neutrons in the whole angle range, i.e., from the angle = 
to = TT in Eq. (^|). Among those neutrons, the ones from the angle ^ = § are most probable because of azimuthal 
symmetry, that is, the number of the incoming neutrons is proportional to 2TTsm9d9. Consequently, we estimate the 
width of the spectrum of the emitted photons from Eq. (Rll) as follows 
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Auji 




y^0.5^AK^ + 12.22Acos6'2 eV 

21.1 eV, for A cos = \/3 and Aii' = 0.005 eV, 



(18) 



where we expanded it at Kq = 0.056 cV and cos^o = 0. For the initial neutron kinetic energy of 0.025 eV, the hne 
width is expected to be 14.0 eV. The above numbers are surprisingly large compared to the error from the crystal 
spectrometer. Here one can expect another line width coming from the uncertainty of the kinetic energy of the 
incident neutrons, expressed by AK in Eq. (p^). But it is very small enough to be neglected as calculated above 
and pointed out in Ref. ||l|. Thus, from the above relation, the broadening of the 7 ray is mainly due to the Doppler 
effect. However most of Ge detectors or the other 7 ray detectors may not resolve the broadened spectrum, so the 
crystal spectrometer is used for precise measurements by determining the Bragg angle as in Ref. [Q . 

In this experimental situation, the broadened 7 ray causes the uncertainties in the measurement of the Bragg angle 
and also in the spectrums of the 7 ray detector which determines the peak positions. The authors of Ref. ||l| have 
measured the 52 Bragg angles between the two centers of the peaks of the broadened spectrums. The small error of 
2.3 eV, which comes from the statistical error in the data of the 52 Bragg angles and the lattice-spacing of crystal, 
is much less than the calculated line width of 14.0 - 21.1 eV. The reason of such a small error may be conjectured 
that the standard error is obtained from the standard deviation by dividing it by the square root of the number 
of measurements and the errors in the center positions of the broadened spectrums are neglected. So the standard 
deviation, which is simply calculated as 2.3 x \/52 — 16.6 eV, is compatible with our calculated line width. Therefore 
if we reduce the line width of the 7 ray by selecting the special angle ^ = § , the measurement of the Bragg angle has 
much smaller error and is much more convincible though it is difficult to reduce the uncertainty of the 7 ray detector 
itself. Actually, the standard deviation of the measured Bragg angle 1.4 x 10~^ rad in Ref. ^ is larger than the 
sensitivity and accuracy of ~ 10~^ arcsec (< 10^^ rad) of the Michelson interferometer, so that it can be reduced up 
to the sensitivity. 

Table | shows the update values of the deuteron binding energy. Among them, only two articles |^,^ have shown 
their experimental geometry. So, we discuss these experimental geometry on a possible Doppler effect. Since these 
experiments seems to be carried out at room temperature, we assume that the kinetic energy of the thermal neutrons 
is 0.025 eV because there are no more informations on the kinetic energy of the incident neutron. In the experiment of 
Ref. 1^ , the incident neutrons seem to be well-thermalized and the neutron source is located to the side of the paraffin 
sample. The neutrons moving perpendicular to the direction of the detected photons are captured dominantly as the 
previous consideration, and the width of the detected photons is calculated as 14.0 eV with using Eq. (|l|). If the 
paraffin sample is more thin and has larger diameter than that of this experiment and exposed to the thermalized 
neutrons moving perpendicular to the detected photon by shielding the sample except the round surface of the sample 
cylinder, one can reduce A cos 9 and the Doppler broadening of the spectrum of the photon. 

In the experiment of Ref. the situation is so complicated to treat incident neutrons as well-thermalized ones. 
There is a possibility that many fast neutrons can be captured through thermalization in the paraffin sample, because 
the source emits much higher energy neutrons than the thermal neutrons. Even though it is difficult to know what 
energy of the neutrons contribute to the capture process dominantly, it is possible to compare this experimental result 
with others in the point of view of the Doppler shift. In Refs. |^,^, the authors have modified the binding energy of 
Ref. [|j on the basis of the same standards of the calibration energies of Refs. as 2.224 628 (15) MeV, and have 

compared it with their measured ones. The corrected value is shifted by 53 eV with respect to the value of Ref. ||] 
and 64 eV with respect to the value of Ref. These shifts occur, when the kinetic energy of the dominant incident 
neutron is more than at least 1.104 eV under the assumption that the incident neutrons are parallel to the direction 
of the detected photons according to our result. The authors of Ref. |^ have also revised their value as 2.224 568 (8) 
MeV on the basis of the standards of Refs. |p|,p| in Ref. In the point of view of our result, this value can not be 
understood easily, on the contrary, this discrepancy is attributed to the problem of the calibration using a common 
system of standard energies. Table ^ shows other similar experiments, in which it is more important to consider the 
Doppler effect in the case of a triton. 

The direction of motion of the thermal neutrons is distributed at random in the previous experiments [Q-Q, and 
Doppler broadening contributes in full to the observed line width. Selecting a certain angle of the incident neutron 
with respect to the emitted photon will lead to a reduced line width and should therefore make it possible to obtain 
even more precise results. Such an experiment would be very challenging. A possible geometry could make use of a 
target sample near the core of a reactor with appropriate shielding from thermal neutrons. Using the difference of the 
penetration lengths for neutrons and photons, shielding the target should not cause much loss in flux of the incident 
neutron beam. As shown by the detailed numbers, it is indispensable to reduce the Doppler broadening for a more 
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precise measurement of the gamma-ray wavelength and the determination of the deuteron binding energy. Careful 
geometric considerations can also explain the discrepancies among the values of the deuteron binding energy reported 
in Refs. 0-0], if calibration problems are properly accounted for. 
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Reference 


Binding energy (oV) 


Greene et al. Q 
Van der Leun et al. 
Vylov et al. [;| 
Vylov et al. \i\ 
Greenwood et al. [^ 
Adam et al 
Michael et alJ^ 


2 224 589.0 (2.2) 
2 224 575 (9) 

2 224 572 (40) 
2 224 568 (8) 

2 224 564 (17) 
2 224 574 (9) 

2 224 579 (13) 


TABLE I. The recently reported values of the deuteron binding energy. The deviations among the values show up at 10 eV 
order, which is just the same order as the Doppler effects contribute (see Eq. (^)). However notice that some of them have 
the error of several electron volts order, which is smaller than the Doppler effects. 


(n,7) reaction product nucleus (MeV) 


7-ray energy (MeV) 


Binding energy (MeV) 


Doppler shift 


13 c 

14 Q 

15 N 


6.250 316 
4.949 319 
8.173 922 
10.829 101 


6.257 268 (24) 
4. 946 329 (24) 
8.176 483 (40) 
10.833 297 (38) 


15.2 cV cose 
2.8 eV cos 6* 

4.2 eV cos e 

5.3 eV cos 9 



TABLE II. These are other nucleus binding energies from (n,7) reactions in Ref. [5] and the possible Doppler shifts according 
to the angle of the incident neutrons with respect to the direction of the detected photons, where we assume that the kinetic 
energy of thermal neutrons is 0.025 eV. 
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